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Abstract 



I review the use of the Heavy Quark Effective Theory in the computation 
of total rates and differential kinematic distributions in inclusive semileptonic 
and radiative B decays. Particular attention is paid to strategies for the 
extraction of Vcb and Vub- 
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I. THE HEAVY QUARK EXPANSION FOR INCLUSIVE B DECAYS 



Why is it important to have a good theoretical understanding of inclusive weak B decays? 
One might be tempted to take the point of view that since the short distance structure of 
these decays is presumably known, all that is left to do is to sort out some messy but 
ultimately unenlightening details of the QCD dynamics. However, by doing so one would 
forget the important fact that while the structure of these decays is indeed well understood, 
the strengths of the weak decay couplings of the b quark are not yet known with sufficient 
precision. The weak decay of the b is governed by two fundamental constants which appear 
as coefficients of currents found in the weak interaction Lagrangian. The larger of the two 
is Vcb, which is multiplies the current C7'^(l — 7^)6; the smaller is Vub, which multiplies 
M7^(l — 7^)6. Respectively, they govern the rates for the inclusive semileptonic decays 
B — » Xcijy and B — > X^^z/. For reasons of theoretical simplicity, I will confine myself to 
the discussion of semileptonic, rather than nonleptonic, weak B decays.0 In addition, I will 
discuss radiative decays of the form B —>■ Xs^iij. These transitions, which are mediated in 
the Standard Model by one-loop penguin operators, may also be studied with the techniques 
to be described here. 

The theoretical tool which is used to analyse semileptonic B decays is the Operator 
Product Expansion (OPE) [0,0. The OPE exploits a fact known as parton-hadron duality, 
which is, roughly speaking, the observation that if the energy release in a decay is large 
{rrif, ^ Aqcd) and the decay is sufficiently inclusive, then it is essentially controlled by 
physics at short distances. Long distance effects, which arise from the B meson bound 
state structure, appear only as subleading corrections. The result is an expansion for the 
semileptonic width F of the generic form 
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where Fq is the decay rate of a free b quark, and Aj are nonperturbative parameters of 
order Aqcd which are independent of in the limit m^, 00. The quantities Aj cannot 
be computed theoretically from first principles (except, perhaps eventually, on the lattice); 
instead, they must be modeled or, preferably, be measured. It is the current state of the art 
to include terms in F up to order l/m\. Those of order 1/mf, and are taken seriously, 

while those of order l/m\ are used (usually crudely) to estimate the residual uncertainties, 
which are typically at the level of a few percent. 

Simultaneously, one may exploit the fact that the the QCD coupling constant is 
perturbative at the scale m^, since as{mh) — 0.2. Since the typical gluon carries a momentum 
set by the scale 3> Aqcd, the theory is in the region where asymptotic freedom applies. 
As a result, each term in the expansion (|I]) is actually a power series in ^^(mft), leading to 
the refined (schematic) form 



^Actually, the analysis of nonleptonic decays is on a less rigorous theoretical footing than that of 
semileptonic decays [j^. 
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It is the present state of the art to compute Ci to order and Cq to "order" a^/3o, where 
/3o = 9 is the first coefficient in the QCD beta function. Such terms, while not formally 
dominant in any nearby limit, are often the largest contribution to the two loop correction [Q. 
Here the largest uncertainties are from the choice of renormalization scale /i, and from 
infrared renormalons which appear, in principle, at high orders in the perturbative expansion. 
The renormalons will be discussed further below. 

To the order which they are computed, there are four input parameters which appear in 
the expansion. First, the QCD perturbation series depends on as{mb), which we will treat 
as known, since the uncertainty associated with its value is comparatively small. Second, 
there are parameters which arise in the OPE. The most important of these is A, defined by 

pole 
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The constituent quark model, along with QCD sum rule estimates, lead one to expect 
200 MeV < A < 700 MeV. At the next order, there appear two matrix elements of dimension 
five operators 

Ai = {B\b{iDfb\B)/2MB, 

\2 = {B\b (^-^a'-'^G,,^ b \B)/6Mb . (4) 

The parameter Ai is related to the negative of the b quark kinetic energy in the B meson; 
models and theoretical prejudice would indicate that it lies in the range < — Ai < 1 GeV^. 
The parameter A2 is related to the matrix element of the leading operator which violates 
heavy spin symmetry, and is proportional to the B — B* mass splitting. Neglecting a small 
radiative correction, A2 ~ 0.12 GeV^. In what follows, we will treat as and A2 as known, 
and A and Ai as parameters which must be determined somehow from experiment. 

In fact, there is an additional subtlety associated with the definition of A, or more pre- 
cisely, the definition of w^^"'^. The problem is that wi^"'*^ is not well defined in QCD. Of 
course, nonperturbatively there is no pole in the fully dressed quark propagator, hence no 
unambiguous definition of m^°^'^. But even within QCD perturbation theory, the asymptotic 
nature of the expansion leads to an ambiguity in the perturbative definition of ml"^'^. In 
particular, one can attempt to sum the "BLM-enhanced" terms proportional to a"/3Q~^, 
where /3o is the first coefficient in the QCD beta function, using Borel resummation tech- 
niques. However, one finds an obstruction in the Borel plane due to an infrared renormalon 
ambiguity, which must be resolved by the choice of a resummation scheme. What this 
means, essentially, is that ml°^'^ has a scheme dependence in its definition, quite analogous 
to the renormalization scale dependence of quantities defined in dimensional regularization 
with MS subtraction. As with the renormalization group, this scheme dependence cancels 
in physical quantities, if the calculation is organized self-consistently. The problem is not 
with the predictivity of the theory, but with the physical interpretation of ml°^^. Since the 
renormalon ambiguity is of order 100 MeV , there is no "preferred" value for m^"^*^ 

within a smaller precision than this. The most important practical point is that if the renor- 
malon ambiguity is not canceled consistently, then there is an irreducible uncertainty in the 
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input A, which it inherits from ra^"^'^ . One of the issues which I will address is the consistent 
treatment of this effect. 



II. INCLUSIVE SEMILEPTONIC B DECAY 

I turn first to the extraction of |V"cfe| from the inclusive decay B Xciv. The expression 
for the semileptonic width is written most conveniently as a function of the ratio of spin 
averaged masses Md/Mb, where Md = {Md + ^Md*)/4. and = (Mb + 3Mb.)/4. These 
physical masses are, in turn, functions of the quark masses mb and iric and the hadronic 
parameters A, Ai and A2, according to 

T\/r pole I A + ^-^2 , 

Mb = ml + A h . . . . 

Zrrify 

MB^=mf' + A-^^^ + .... (5) 

zruh 



The result for the semileptonic width is P,^,|9|JT0| 



T{B XJu) = ^^i^^ / ( Mo/Mb) |i - 1.54^ - (1.43/3o + c) (^^' 



1927r'^ ^ -^1 TT V TT 

A / as\ A^ Ai A2 

>— 1 - 0.87— - 0.95—^ - 3.18-^ + 0.02-^ 

Mb\ n J Ml Ml Ml 



where /{Md/Mb) = 0.369 is a phase space factor, and the two loop coefficient c is known to 
be small . An important feature of the QCD expansion is that there is also a renormalon 
ambiguity in the perturbation series in Eq. @; that is, there is a scheme dependence in 
the summation of the series of terms proportional to a^P^'^. However, one can show 
that this ambiguity precisely compensates the renormalon ambiguity in A, so the scheme 
dependence cancels in the physical width The mechanism is similar to the cancelation 
of renormalization scale dependence in physical quantities. 

Once T{B ^ Xclu) has been measured, the expression (||) can be used to extract \Vcb\- 
There is a variety of sources of uncertainty in this determination. First, both the perturbative 
QCD expansion and the operator product expansion are truncated, and we must estimate 
the size of the omitted terms. For the radiative corrections, the leading unknown term is of 
order (as/vr)^, and is expected to be no larger than a few percent. For the nonperturbative 
corrections, the next terms are of order (A/Mb)^, also at the level of a percent or so. Second, 
there is uncertainty in the values of the parameters to be inserted in the terms which have 
been calculated. By far the most important of these is the uncertainty in A, or equivalently 
in ml"^^, because it enters the expansion already at first order. For example, an uncertainty 
in A of 200 MeV implies an uncertainty in \ Vcb\ of approximately 10%. This is the dominant 
source of error in the determination of \Vcb\- 

The cleanest and best way to reduce this uncertainty is to determine A directly from the 
semileptonic decays themselves. This is possible because one can compute not only total 
decay rates with the OPE, but kinematic distributions as well. While there exist interest- 
ing pathologies in certain of these distributions, low moments of kinematic observables are 
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almost always well behaved. In particular, it is useful to consider the following two types of 
quantities: 

(i) In the decay B Xciiy, let sh be the invariant mass of the hadronic state X^. Then 
the moments of the form {{sh — ^ d)"") have an OPE and a perturbative QCD expansion. 
Constructed in this way, the first moment starts at order AMb, while the second starts 
at order A^. Since terms of the order are sources of uncertainty, the first moment is 
substantially more reliable theoretically than is the second. These moments have been 
studied in detail in Refs. |1^,|T3[. 

(m) One may also study the differential distribution dT/dEi and its moments {E^). 
Again, these quantities have an OPE and a perturbative QCD expansion. They have been 



studied in detail in Refs. |T^ , [T5[| 



Using a calculation analogous to that of the total semileptonic decay rate, these kine- 
matic moments have expansions in terms of a^. A, Ai and A2. As with the total rate, 
the renormalon ambiguity in the perturbative series cancels that in A. Thus, if one writes 
r(S — > Xciiy) in terms of one of these moments, eliminating A, the renormalon ambiguity 

cancels automatically in the resulting expression. For example, if we eliminate A in favor of 
2 

Si = {sh — M^) and solve for \Vcb\, we find 
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+ 0.68-^ 1 + 0.70— + 0.38-^ - 2.00-^ + 1.14-^ + . . . . (7) 

M| V 71 J M% Ml Ml ] ^ ^ 

In this expression, the two loop contributions and the uncertainty due to Ai are both at the 
level of one percent. The effect of the leading infrared renormalon has been canceled, and 
one expects the remaining perturbation series in to be much better behaved than in the 
original expansion (^.0 Of course, one then needs an accurate measurement of Si or another 
kinematic moment. In Ref. it is shown that the theoretical sensitivity of Si and {Eg) 
to A are roughly the same, once corrections of order {K/MbY ^^nd necessary kinematic cuts 
are included. Either way, such an analysis is the strategy which will be required to reduce 
further the experimental error of \Vcb\- 



III. CHARMLESS INCLUSIVE SEMILEPTONIC B DECAY 

I turn now to charmless semileptonic B decay, and the inclusive width for B X^iv- 
The partial decay rate YiB ^ Xulv) is calculable as before, but unfortunately it is not 
very useful phenomenologically. The reason is that |Ktf,|/|Vcfe| ~ 0.1, so approximately 99% 
of semileptonic B decays have charm in the final state. The only way to reject this huge 
background is by making kinematic cuts which exclude charmed final states unambiguously. 

The simplest method experimentally, which has been used extensively in the past, is 
to require that the charged lepton in the decay have an energy larger than approximately 



^An explicit example of how this improvement works, up to order a^, is given in Ref. ||T 
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2.2 GeV, beyond the kinematic endpoint for decays to charm. This constitutes only a small 
tail of the spectrum dT/dEi. In order to extrapolate to the complete partial width, it 
is necessary to have an accurate theoretical understanding of the shape of dT/dEi in the 
endpoint region. Unfortunately, while the OPE can be used to predict dT/dEi over most 
of the available phase space, it is actually an expansion in powers of A/ {Mb — 2E£) rather 
than A/Mb M^M- In the endpoint region, where Ee ~ Mb/2, the OPE fails to converge 
reasonably and cannot be used for an accurate analysis. 

As a result, the extraction of \Vub\ by this method is extremely model dependent. Two 
approaches are used to address the issue. The first is inclusive, in which one performs the 
calculation as before but replaces the higher order terms in the OPE with a "shape function" , 
which is an ansatz which resums formally an infinite set of terms in the series. This shape 
function is usually modeled by giving the b quark a Gaussian momentum distribution ||17|| . 
The second is to treat the endpoint region as a sum over decays to individual exclusive 
final states, such as Triu, piu, which are then analysed with a constituent quark potential 



model 18,19 



Each of these approaches has obvious deficiencies and misses important physics. Fur- 
thermore, even within each model there are uncertainties which cannot be quantified or even 
estimated, since these models are in no sense systematic approximations to QCD. There is 
an important moral to be drawn. In a situation such as this, analysing additional models 
does not lead to better accuracy in the determination of \Vub\- Nor can the uncertainty in 
\Vub\ due to the use of models be estimated sensibly simply by surveying the models currently 
available on the market. The model dependence in \Vub\ is impossible to quantify reliably, 
and it continues to be underestimated consistently throughout the literature on this subject. 

However, with the development of techniques for neutrino reconstruction, a new type 
of analysis has become possible. If the energy and momentum of the missing u can be 
reconstructed kinematically, then it is possible to measure inclusively the invariant mass sh 
of the final state Xu- Charmed final states then can be rejected by requiring < < M^,, 
where is an experimental cut which may have to be less than M^ because of details of 
the neutrino reconstruction procedure. The advantage of this approach is that most of the 
final states in B X^ii^ will pass such a cut, if A^ is reasonably close to M^. One still 
must correct for the small "leakage" of rate to larger values of sh, but, in contrast to the 
previous analysis, one is not extrapolating most of the rate from a small tail. Hence the 
larger A^ can be made to be, the smaller the effect of this tail on the measurement of \Vub\, 
even if the tail is not modeled particularly well. 

A careful analysis |2^ (see also Ref. [^) reveals that the fraction of rate with sh > A^ 
depends crucially on the value of A. This is true for two reasons. First, at tree level in 
the parton model, the spectrum dT/dsn cuts off at sh = Amb, considerably below M|,. 
Thus, gluon bremmstrahlung plays an important role in creating states of large sh, and 
this bremmstrahlung has a Sudakov-type double logarithmic singularity at sh = Amb- The 
position of this (integrable) singularity controls the strength of the tail of the spectrum. 
Second, the nonperturbative corrections from the OPE have the effect of smearing the parton 
level spectrum by an amount of order Arub- When this smearing is modeled in the ACCMM 
model, it turns out that the dominant sensitivity of the size of the tail is to the value of A; 
once A is fixed, the dependence of the leakage on other parameters, such as Ai, is minimal. 

Hence it is important, once again, to measure A using one of the inclusive methods 
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discussed in the previous section. An accurate determination of this parameter is crucial for 
the accurate determination of \Vub\- In addition, we should hope that the value of A turns 
out to be rather small, and that the experimenters can make close to M|,, so that the 
separation between Am^ and will be as large as possible. From the analysis of Ref. ||20|| , 
one may estimate the eventual theoretical accuracy under a variety of hypothetical 

scenarios: (i) if A^ ~ M|, and A = 400 ± 100 MeV, then ~ 10%; (n) if A^ ~ 1.5 GeV^ 

and A = 400 ±100 MeV, then S\Vub\ ~ 30%; (Hi) if A^ - 1.5 GeV^ and A = 200 ±100 MeV, 
then 5\Vub\ ~ 10%; and (iv) if A^ ~ and A = 600 ± 100 MeV, then 5\Vub\ ~ 50%. With 
luck, this method eventually will yield the most accurate value of \Vub\ available, a value 
largely free of model dependence. 



IV. RADIATIVE B DECAY 

Finally, I turn to the extraction of information about short distance dynamics from the 
inclusive radiative decay B Xs'y. Away from the region in which M{Xs) ~ M{J/ip), this 
decay is mediated primarily by a transition magnetic moment operator of the form sa^^Ffj^^b. 
In the Standard Model, this operator is induced by one loop GIM-violating effects, and its 
coefficient is small. The total branching fraction for this decay in the Standard Model is at 
the level of 10~^. One the other hand, the coefficient of this operator, and hence the decay 
rate in this channel, is substantially larger in many models of new physics, especially those 
that postulate new flavor dynamics at the TeV scale. In order for an observed enhancement 
of 5 — *■ Xs'^ to serve as a signal of new physics, we must know how to predict the physical 
inclusive rate from a known quark operator such as sa'^'^F^yb. 

In fact, the situation is quite analogous to that for B X„£z/, in that the total rate for 
B Xs'~f can be computed reliably with the OPE but is not quite useful phenomenologi- 
cally ||22|,^ . Once again, the difficulty is the much larger number of decays to charmed final 



states. The most straightforward way to search for B Xs'j is to look for the relatively 
hard 7, but there is an important background from the process B Dn^ D'-f-f, where 
one of the photons is missed. If one looks at the photon energy spectrum dT/dE^, it is 
unambiguously from the final state X^'f only for > 2.2 GeV. As a result, one is again in 
the position of extrapolating the bulk of the spectrum from a kinematic tail. For reasons 
similar to those which plague B — > the OPE is not convergent in the region to which 

one is restricted. 

However, in some respects the situation is not quite as severe as for B Xuiu. At 
tree level in the parton model, the photon is monochromatic with E^ ~ 2.4 GeV, above 
the kinematic cut. Radiative corrections and nonperturbative bound state effects smear the 
spectrum about this point, but most of the decay rate remains at large values of E^. In 
particular, even a relatively small downward shift in the value of Ej at which one must cut 
can reduce substantially the fraction of the rate which must be inferred by extrapolation. 
At present, however, the state of the art is still to use an inclusive model in this analysis, 
based on an Gaussian ansatz for the initial momentum distribution of the b quark in the B 
meson 



24 



What is most important for the future is to aim for a measurement of dV/dE^ down to 
values of E^ as small as possible. Also, an accurate measurement of the shape of the spectrum 
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in the region > 2.2 GeV would help substantially to reduce the model dependence of 
the result, by providing information to constrain the models which are used. An accurate 
measurement of the spectrum in this region also could provide information on the shape 
function which governs the endpoint of the lepton energy spectrum in B ^ X^C-V [|l6| , p5| , or 
could be used in an extraction of A and Ai which is complementary to that from the decay 

B xju M,m. 



V. SUMMARY 



Inclusive B decays are indispensable for the accurate extraction of \Vcb\, \Vub\, and the 
coefficient of the operator responsible for the decay b — *• sj. While the theoretical treatment 
of such decays has developed into a mature field over the past five years, a number of 
subtleties remain. For the extraction of \Vcb\ from the dominant decay B —>■ Xcii', the 
most important theoretical uncertainties are those in the inputs A (or rrib) and (to a lesser 
extent) Ai, and from the truncation of the perturbative expansion at order a^. However, 
the infrared renormalon relates an ambiguity in the definition of A to an ambiguity in the 
resummation of the perturbative series, in such a way that the consistent elimination of A 
from physical quantities also helps to control the bad behavior of the perturbation series. 
Explicit strategies for extracting A from experiment are discussed. 

For the rare decays B —>■ Xu^v and B — >■ Xs7, more serious theoretical uncertainties are 
induced by the severe kinematic cuts which must be imposed to eliminate the background 
from B —>■ Xciiy. A strategy for extracting \Vub\ from B Xu^v-, which relies on the 
experimental technique of neutrino reconstruction, has been developed. Future improvement 
in the theoretical understanding of inclusive semileptonic and rare B decays will depend on 
{i) the direct measurement of A and Ai in inclusive decays themselves, a program which is 
already underway; and {ii) application of the loosest kinematic cuts possible in rare decays, 
which will help reduce the model dependence of the results. In the short and in the long 
term, the future in this field is bright; ultimately, the prospects are excellent for obtaining 
accurate and reliable insight into the CKM matrix and other short distance physics from 
inclusive semileptonic and radiative B decays. 
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